Toroidal dipole resonances in the relativistic random phase 

approximation 

D. Vretenar^'^, N. Paar^, P. Ring^, and T. Niksic^-^ 

1. Physics Department, Faculty of Science, University of Zagreb, 10000 Zagreb, Croatia 

2. Physics Department, TU Munich, D- 857 4-8 G arching, Germany 

O ; (February 8, 2008) 

(N 

C^ ' The isoscalar toroidal dipole strength distributions in spherical nuclei are 

OO \ calculated in the framework of a fully consistent relativistic random phase 

^ _ approximation. It is suggested that the recently observed "low-lying compo- 

T^ ' nent of the isoscalar dipole mode" might in fact correspond to the toroidal 
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^D ' giant dipole resonance. Although predicted by several theoretical models, the 
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existence of toroidal resonances has not yet been confirmed in experiment. 
The strong mixing between the toroidal resonance and the dipole compres- 
sion mode might help to explain the large discrepancy between theory and 
experiment on the position of isoscalar giant dipole resonances. 
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j^ ■ In addition to data on isoscalar giant monopole resonances (IS GMR) in spherical nuclei, 

it is expected that information on nuclear matter incompressibility should also be obtained 
in studies of the isoscalar dipole mode. The isoscalar giant dipole resonance (IS GDR) is a 
second order effect, built on "ihuo, or higher configurations. It corresponds to a compression 
wave traveling back and forth through the nucleus along a definite direction. Recent data 
on IS GDR obtained by using inelastic scattering of a particles on ^osp]^ ^^ o^j^j qj^ ^^Zr, 
^^^Sn, ^^^Sm, and ^osp]^ gj^ have been analyzed in the non-relativistic Hartree-Fock plus 
RPA framework P], and with relativistic mean- field plus RPA (RRPA) calculations [Q. Both 
analyses have shown that: (a) there is a strong disagreement between theory and the reported 
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experimental data on the position of the IS GDR centroid energies, and (b) calculations 
predict the splitting of the IS GDR strength distribution into two broad structures, one in 
the high-energy region above 20 MeV, and one in the low-energy window between 8 MeV 
and 14 MeV. Effective interactions, both non-relativistic and relativistic, which reproduce 
the experimental excitation energies of the IS GMR, predict centroid energies of the IS 
GDR that are 4 — 5 MeV higher than those extracted from small angle a-scattering spectra. 
This disagreement between theory and experiment is an order of magnitude larger than for 
other giant resonances. Another puzzling result is the theoretical prediction of a substantial 
amount of isoscalar dipole strength in the 8 — 14 MeV region. In Ref. |^ we have shown 
that the RRPA peaks in this region do not correspond to a compression mode, but rather 
to a kind of toroidal motion with dynamics determined by surface effects. In a very recent 
article ^, Clark et al. reported new experimental data on the isoscalar dipole strength 
functions in ^°Zr, ^^^Sn, and ^°^Pb, measured with inelastic scattering of a particles at small 
angles. They found that the isoscalar El strength distribution in each nucleus consists of a 
broad component at E^ ^ llA/A^^^ MeV containing approximately 100% of the El EWSR, 
and a narrower one at Ex ~ 72/^4^/^ MeV containing 15 — 28 % of the total isoscalar El 
strength. The higher component is identified as the El compression mode, whereas the lower 
component may be the new mode predicted by the RRPA analysis of Ref. [^]. In the present 
work we suggest that the observed low-lying El strength may correspond to the toroidal 
giant dipole resonance (TGDR). 

The role of toroidal multipole form factors and moments in the physics of electromag- 
netic and weak interactions has been extensively discussed in Refs. and 0. They appear 
in multipole expansions for systems containing convection and induction currents. In par- 
ticular, the multipole expansion of a four-current distribution gives rise to three families of 
multipole moments: charge moments, magnetic moments and electric transverse moments. 
The later are related to the toroidal multipole moments and result from the expansion of 



the transverse electric part of the current. The toroidal dipole moment, in particular, de- 
scribes a system of poloidal currents on a torus. Since the charge density is zero for this 
configuration, and all the turns of the torus have magnetic moments lying in the symmetry 
plane, both the charge and magnetic dipole moments of this configuration are equal to zero. 
The simplest model is an ordinary solenoid bent into a torus. 

Vortex waves in nuclei were analyzed in a hydrodynamic model ^j. By relaxing the 
assumption of irrotational motion, in this pioneering study solenoidal toroidal vibrations 
were predicted, which correspond to the toroidal giant dipole resonance at excitation energy 
E^ ^ (50 — 70)/A^/^ MeV. It was suggested that the vortex excitation modes should appear 
in electron backscattering. In the framework of the time- dependent Hartree-Fock theory, 
the isoscalar 1^ toroidal dipole states were also studied by analyzing the dynamics of the 
moments of the Wigner transform of the density matrix [^ . 

In this Letter the toroidal dipole strength distributions are calculated in the relativistic 
random phase approximation (RRPA). The RRPA represents the small amplitude limit of 
the time-dependent relativistic mean-field theory |1^. A self-consistent calculations ensures 
that the same correlations which define the ground-state properties, also determine the be- 
havior of small deviations from the equilibrium. The same effective Lagrangian generates the 
Dirac-Hartree single-particle spectrum and the residual particle-hole interaction. In a num- 
ber of recent applications [3,|riHl5| it has been shown that, by using effective Lagrangians 
which in the mean-field approximation provide an accurate description of ground-state prop- 
erties, excellent agreement with experimental data is also found for the excitation energies of 
low-lying collective states and of giant resonances. Two points are essential for the successful 
application of the RRPA in the description of dynamical properties of finite nuclei: (i) the 
use of effective Lagrangians with non-linear terms in the meson sector, and (ii) the fully 
consistent treatment of the Dirac sea of negative energy states. In particular, in Ref. [|1^] it 
has been shown that configurations which include negative-energy states have an especially 



pronounced effect on isoscalar excitation modes. 

In Fig. |l](a) we display the RRPA toroidal dipole strength distribution for ^°^Pb: 

R(E) = y B^=^(E1, ir -^ Of) — -, 

where F is the width of the Lorentzian distribution, and 
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For the strength distributions in Fig. |T] the width is F = 1.0 MeV. The isoscalar toroidal 
dipole operator is defined 
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In the relativistic framework the expression for the the isoscalar baryon current reads 
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where the summation is over all occupied states in the Fermi sea. The resulting toroidal 
dipole operator is 
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where Y Wfj, denotes a vector spherical harmonic, and a are the Dirac ct-matrices. The 
calculations have been performed with the self-consistent Dirac-Hartree plus relativistic 
RPA. The effective mean-field Lagrangian contains nonlinear meson self-interaction terms, 
and the configuration space includes both particle-hole pairs, and pairs formed from hole 
states and negative-energy states. The inclusion of the term — < r^ >(, i iqu i^ ^^^ operator 
ensures that the TGDR strength distributions do not contain spurious components that 
correspond to the center-of-mass motion [§]. In Fig. |l|(a) we compare the toroidal strength 
distributions calculated without (dashed) and with (solid) the inclusion of this term in the 
operator. 



The projection of spurious center-of-mass motion components can be also performed by 
subtracting the spurious transition current, following a procedure similar to that adopted in 
Refs. 0,0 for the IS GDR. The toroidal dipole strength distribution can be written as 
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where T{r) is the vector toroidal operator, and Sj{f) is the transition current 
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The radial functions j-(r) and j+(r) are defined in Ref. |T^. We have verified that identical 
strength distributions (solid curve in Fig. |l|) are obtained when: (a) the transition current 
(^) is used and the toroidal operator is corrected by including the term — < r^ >q Y iq^, 
or (b) this term is not included in the operator and the spurious component is subtracted 
from the transition current at each energy 



6j{r) - apoe;. 
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The second term in this expression is the spurious transition current |]18[. po is the ground- 
state density and Cz denotes the unit vector in the direction of the center-of-mass motion. 
The energy dependent coefficient a is determined by the condition that the integral of the 
transition current over the nucleus should vanish at each energy 



d^r (6j{r) — ap^'t 
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The strength distributions in Fig. ^ have been calculated with the NL3 [|T^ effective in- 
teraction. In Ref. [^ it has been shown that isoscalar giant monopole resonances calculated 
with this effective force (-ft'nm = 271.8 MeV) are in excellent agreement with experimental 
data, and in Ref. 0] this interaction was used in the RRPA analysis of the IS GDR. By 
using effective interactions with different values of the nuclear matter compressibility mod- 
ulus, it was shown that only the high-energy (above 20 MeV) portion of the isoscalar dipole 



strength distribution corresponds to a compression mode. The same effect is observed for 
the toroidal strength function: the position of the peaks in the low-energy region (below 20 
MeV) depends only weakly on the incompressibility, while the structure in the high-energy 
region is strongly affected by the choice of the compression modulus of the interaction. In 
Fig. |T](b) we plot the strength function of the isoscalar dipole compression operator ^ 

Ql;' = E 70 (r^ - ^ < r^ >o r) ¥,,{6,, ^,). (10) 

j=i -J 

We note that both dipole strength distributions, toroidal in Fig. |l](a) and compressional in 

Fig. |l](b), display two broad structures: one at low energies between 8 and 15 MeV, and 

one in the high-energy region 25 — 30 MeV. Obviously, one could expect a strong coupling 

between the two isoscalar 1~ modes. This coupling becomes even more evident if one rewrites 

the expression in square brackets of the toroidal operator (^ as [§] 

Vx (fx V)(r=^-- <r2>or) Fi^, (11) 

and compares it with the isoscalar dipole operator of the compression mode (0). The 
relative position of the two resonance structures will, therefore, depend on the interaction 
between the toroidal and compression modes. 

In Fig. we display the RRPA toroidal dipole strength distributions in ^°Zr, ^^^Sn, and 
^°^Pb, calculated with the NL3 effective interaction. In all three nuclei a broad, strongly 
fragmented structure is found in the low-energy region, where "the low- lying component of 
the IS GDR" has been observed 0. The toroidal strength distributions in this region should 
be compared with the experimental centroid energies of the "low-lying component" [Q: 
16.2±0.8 MeV for ^^Zr, 14.7±0.5 MeV for ^^^Sn, and 12.2±0.6 MeV for ^o^Pb. The calculated 
peaks in the high-energy region, on the other hand, correspond to the compression mode. 
The dynamics of the solenoidal toroidal vibrations is illustrated in Fig. ^, where we plot the 
velocity fields for the four most pronounced peaks of the toroidal dipole strength distributions 
in ^^^Sn (see Fig. ^. The velocity distributions are derived from the corresponding transition 
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currents, following the procedure described in Ref . |T^ . A vector of unit length is assigned 
to the largest velocity. All the other velocity vectors are normalized accordingly. Since the 
collective flow is axially symmetric, we plot the velocity field in cylindrical coordinates. The 
2;-axis corresponds to the symmetry axis of a torus. We note that the two lowest peaks at 
8.82 MeV and 10.47 Mev are completely dominated by vortex collective motion. The velocity 
fields in the {z, r±) plane correspond to poloidal currents on a torus with vanishing inner 
radius. The poloidal currents determine the dynamical toroidal moment. The high-energy 
peak at 30.97 MeV displays the dynamics of dipole compression mode. The "squeezing" 
compression mode is identified by the flow lines which concentrate in the two "poles" on 
the symmetry axis. The velocity field corresponds to a density distribution which is being 
compressed in the lower half plane, and expands in the upper half plane. The centers of 
compression and expansion are located on the symmetry axis, at approximately half the 
distance between the center and the surface of the nucleus. Finally, the intermediate peak 
at 17.11 MeV clearly displays the coupling between the toroidal and compression dipole 
modes. A very similar behavior of the velocity distributions as function of excitation energy 
is also observed for ^°Zr and ^''^Pb. 

We suggest, therefore, that the recently observed "low-lying component of the isoscalar 
dipole mode" might in fact correspond to the toroidal giant dipole resonance. By em- 
ploying the fully consistent relativistic random phase approximation, in Ref. [Q and in the 
present analysis we have shown that the toroidal dipole strength is concentrated in the low- 
energy region around 10 MeV, while the isoscalar dipole excitations in the high-energy region 
above 20 MeV correspond to the "squeezing" compression mode. States in the intermediate 
region display a strong mixing between the two dipole resonances. The pronounced coupling 
between the toroidal resonance and the IS GDR, predicted by the RRPA calculations, might 
also explain the strong discrepancy between theory and the experimental position of the IS- 
GDR centroid energies [l5|,p!6|,|3|J^ . Namely, the interaction causes a repulsion of the two 1~ 



isoscalar modes, i.e. it pushes to higher energies the IS GDR and pulls to lower energies 
the TGDR. This effect would explain the observation of Ref. ^ that: "The centroids of the 
higher (compression) mode calculated with interactions which reproduce GMR energies are 
about 4 MeV higher than the experimental centroids while the calculated centroids for the 
lower mode lie 1-2 MeV below the experimental values." 
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FIG. 1. (a)Toroidal dipole strength distributions in ^'^'^Pb, calculated without (dashed) and 
with (solid) projection of spurious center-of-mass components, (b) IS GDR strength distributions 
in 208pb. 

FIG. 2. Toroidal dipole strength distributions in ^^Zr, ^^^Sn and ^^^Pb, calculated with the NL3 
effective interaction. 

FIG. 3. Velocity distributions for the most pronounced dipole peaks in ^^^Sn (see Fig. y). The 
velocity fields correspond to the peaks at 8.82 MeV (a), 10.47 MeV (b), 17.11 MeV (c), and 30.97 
MeV (d). 
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